On  a  Diophantine  Equation  of  Stroeker 

F.  Luca,  P.  Stanica,  and  A.  Togbe 
October  27,  2008 


Abstract 

In  this  paper,  we  prove  that  there  are  infinitely  many  positive 
integers  N  such  that  the  Diophantine  equation  ( x 2  +  y)(x  +  y2)  = 
N(x  —  y)3  has  no  nontrivial  integer  solution  (x,y). 


1  Introduction 

Let  N  be  a  nonzero  integer.  In  [6],  Stroeker  investigated  the  Diophantine 
equation 


{x2  +  y)(x  +  y2)  =  N(x  -  y)3,  (x,y)  e  Z2.  (1) 

It  clearly  suffices  to  consider  the  case  when  N  >  0  since  if  (x,  y )  is  a  solution 
of  equation  (1)  above,  then  (y,x)  is  a  solution  of  (1)  with  N  replaced  by 
—N.  We  shall  only  consider  solutions  (x,y)  such  that  xy  ^  0.  Following 
Stroeker,  we  refer  to  such  solutions  as  proper.  Note  that  equation  (1)  always 
admits  the  solution  x  —  y  —  —  1.  This  will  be  referred  to  as  the  trivial  so¬ 
lution.  Stroeker  proved  that  if  (x,y)  is  any  proper  solution  to  equation  (1), 
then  max{|x|,  |y|}  <  N 3  if  N  >  4.  He  also  showed  that  if  N  >  1  is  odd, 
then  equation  (1)  has  at  least  one  non-trivial  proper  solution  and  if  addi¬ 
tionally  27 N2  —  2  is  a  square,  then  equation  (1)  has  at  least  5  non-trivial 
proper  solutions.  In  [4],  Mqkowski  pointed  out  a  connection  between  (1) 
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and  Fibonacci  numbers.  In  fact,  he  showed  that  if  N  —  Fk_2Fk_i1  then 
(x,y)  =  {FkFk+l)-Fl+1)  (for  k  even)  and  (. x,y )  =  (F%+1,—FkFk+i)  (for  k 
odd)  are  solutions  of  (1).  Moreover,  he  asked  the  following  question:  ”  Do 
there  exist  infinitely  many  positive  (even)  integers  N  such  that  equation  (1) 
has  only  the  trivial  solution”  ?  Computations  in  the  range  1  <  N  <  51 
showed  that  if  N  e  {8, 10, 12, 14, 16,  20,  24,  26,  28,  30,  36, 44, 48},  then  equa¬ 
tion  (1)  has  no  proper  non-trivial  solution. 

Our  main  result  is  an  affirmative  answer  to  this  question. 

Theorem  1.  There  are  infinitely  many  positive  integers  N  such  that  equa¬ 
tion  (1)  has  no  non-trivial  proper  solution  (x,y). 

Throughout  this  paper,  we  use  the  Vinogradov  symbols  3>  and  -C  as  well 
as  the  Landau  symbol  O  with  their  usual  meanings.  The  constants  implied 
by  them  are  absolute. 

2  Proof  of  Theorem  1 

We  let  X  be  a  large  positive  real  number.  We  put  N  =  2 p,  where  p  G 
(X/2,X)  is  a  prime.  Stroeker  showed  that  for  any  non-trivial  solution  (x,y) 
of  equation  (1)  there  are  integers  u  >  2,  v  >  1  and  £  0  such  that 

2x  =  v  —  u  +  £  +  1,  2y  =  v  —  u  —  £  +  1,  uv  —  N£, 

and 

(u  +  v  —  £)2  =  4  (u  —  l)(v  +  1)  +  1. 

Stroeker  showed  that  max{«,  r}  <  2N3/3  <  6A"3  if  p  >  5,  and  showed  also 
that  both  formulae 

u2{v  -  N)2  -  2 uN(v2  +  vN  +  2 N)  +  N2(v  +  l)(u  +  3)  =  0  (2) 

and 

v2(u  -  N)2  -  2 vN{u2  +  uN  —  2 N)  +  N2{u  -  l)(w  -  3)  =  0  (3) 

hold.  Furthermore,  both  u(v  —  N)2/N  and  v{u  —  N)2/N  are  integers.  If 
v  =  N,  then  4w  =  N  +  3,  which  is  impossible  since  N  +  3  is  odd.  Thus, 
v  N .  The  same  argument  shows  that  u  N .  Since  N  =  2 p,  it  follows 
easily  from  the  fact  that  u{v  —  N)2/N  is  an  integer  that  p  \  uv. 
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Let  us  assume  that  p  \  v  ( Case  1).  Write  v  =  pX,  where  A  fz  2.  Note  that 
A  <  2N3/(3p )  <  6A"2.  Replacing  v  by  A p  in  equation  (2)  and  simplifying  a 
factor  of  p2,  we  get 

u2(  A  —  2)2  —  4u(X2p  +  2A  p  +  4)  +  4(A  p  +  l)(Ap  +  3)  =  0.  (4) 

(If  instead  p  \  u,  then  u  =  pX,  and  equation  (3)  simplified  by  a  factor  of  p2 
becomes  v2(X  —  2)2  —  4v(X2p  +  2Xp  —  4)  +  4(A p  —  1)(A p  —  3)  =  0.  We  shall 
refer  to  this  as  Case  2.) 

The  relation  (4)  of  Case  1  can  be  rewritten  as 

U 2  -  2AR2  =  4  -  2A,  (5) 

where 

U  =  2Xp  +  4  —  uX  —  2  u,  V  —  2u  —  1. 

(If  u  =  pX,  then  the  corresponding  Pell  equation  is  also  U2  —  2XV2  =4  —  2A, 
with  U  =  2Xp  -  4-vX-  2v  and  V  —  2v  +  1.) 

Let  T  :=  T( X)  <  6X2  be  some  parameter  depending  on  X  and  tending  to 
infinity  with  X,  to  be  made  more  precise  later.  Assume  that  A  <  T  is  given. 
If  2A  is  a  square,  then  U  —  V2XV  and  U  +  y/2XV  are  two  divisors  of  4  —  2A 
whose  product  is  4  — 2A.  Hence,  the  pair  (U,  V )  can  be  determined  in  at  most 
2r(|2A— 4|)  ways,  where  for  a  positive  integer  m  we  write  r(m)  for  the  number 
of  its  positive  divisors.  Note  that  the  triple  (A,  U,  V )  determines  p  uniquely. 
Thus,  the  number  of  possibilities  for  the  prime  p  when  2A  is  a  square  is 
•C  t(|2A  —  4|).  Assume  now  that  2 A  is  not  a  square.  It  is  then  well-known 
from  the  theory  of  quadratic  fields  that  there  exist  t  fundamental  positive 
integer  solutions  (Ui,  Vi), . . . ,  (Ut,  Vt)  of  equation  (5)  in  the  following  sense: 
if  (I/,  V)  is  any  positive  integer  solution  of  equation  (5),  then  U  +  \[2XV  = 
(Ui  +  v/2AId)Cm  holds  for  some  i  —  1, . . . ,  t  and  some  nonnegative  integer  m, 
where  we  put  £  for  the  fundamental  unit  of  the  real  quadratic  held  Q[\/2A]. 
Since  £  >  (1  +  Vb)/2  and  max{|t/|,  |V|}  -C  X 5,  it  follows  that  m  <C  log  A". 
It  is  known  that  the  number  t  of  fundamental  solutions  to  equation  (5)  is 
•C  r(|2A  —  4|).  Observe,  as  before,  that  the  triple  (A ,U,V)  determines  p 
uniquely.  Hence,  the  total  number  of  primes  p  that  can  arise  in  this  way 
when  1  <  A  <  T  and  A  2  is 

«  51  r(|2A  —  4|)  log  A  «  T(logT)(log  A)  «  T(log  A)2.  (6) 

1<A<T 

A/2 
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We  now  assume  that  A  G  (T,6X2).  Stroeker  showed,  under  the  condition 
v  7^  N,  which  is  the  case  for  the  positive  integers  N  we  are  considering,  that 
there  exists  a  positive  integer  z  such  that 

0 v 2  +  vN  +  2N)2  -  (v  -  N)2(v2  +  Av  +  3)  =  z2,  (7) 

and 

N(v2  +  vN  +  2N±z) 

7/  =  - - - - 

(v  -  N )2 

If  one  looks  at  (3)  instead,  one  gets,  under  the  assumption  u  ^  N  which  is 
the  case  for  us,  that  there  exists  an  integer  w  such  that 

(u2  +  uN  —  2N )2  —  (u  —  N)2[u2  —  Au  +  3)  =  w 2, 


with 


v  = 


N  ( u 2  +  uN  —  2N  ±  w) 


(■ u  -  N)2 

Returning  to  formula  (7),  let  us  observe  that 

2 


(y2  +  vN  +  2N )  2  =  vA  [l  +  O  (^-  )  )  =,4(l  +  0(^ 


V  v 


X 


while 


TV 

(v  —  N)2{y2  +  Av  +  3)  =  v4  I  1  +  O  I  —  )  )  I  1  +  O  I  - 


=  u4(l  +  0(- 


Hence, 


z2  =  v4(l  +  0[i  =01-  , 


A 


A 


giving 


Thus, 


N(l  +  v/N  +  2N/v2  ±z/v2) 
(1  -  N/v )2 


N  [1  +  0 


N  +  0 


A 


Let  c  be  the  constant  implied  by  the  above  Landau  symbol  and  let  Y  = 
cX/T1/2.  Then  u  =  N  +  m,  where  0  <  \m\  <  Y .  Replacing  in  relation  (3) 
the  variable  u  by  2p  +  m  and  simplifying  a  factor  of  p2,  we  get  the  relation 

A  2m2  —  4A((2p  +  m )2  +  2(2  p  +  m)p  —  4 p)  +  4(2  p  +  m  —  l)(2p  +  m  —  3)  =  0. 

(In  Case  2,  we  replace  v  =  N  +  m  =  2p  +  m  in  (2)  and,  after  simplifying  a 
factor  of  p2,  we  obtain  A 2m2  —  4A((2p  +  m )2  +  2(2p  +  m)p  +  4 p)  +  4(2p  +  m  + 
1)(2 p  +  m  +  3)  =  0.)  Multiplying  both  sides  of  the  above  relation  by  1  —  2A 
and  regrouping  it  we  get 

bb2  —  2A3m2  —  A2(12m  —  4)  +  8A  +  4,  (8) 

where  W  =  4(1  —  2A )p  +  2 (m  —  2)  —  A(3m  —  2).  (In  a  similar  way,  in  Case 
2,  we  obtain  bb2  =  2A 3m2  +  A2 (12m  +  4)  +  8A  +  4,  where  bb  =  4(1  —  2A)p  + 
2 (m  +  2)  —  A (3m  +  2).)  It  is  easy  to  check  that  the  two  variable  polynomial 

P(A,  M)  =  2A3M2  -  A2(12M  -  4)  +  8A  +  4 

is  irreducible  as  a  polynomial  in  C[A,  M].  Indeed, 

AP(A,  M)  =  2((A 2M  -  3A)2  +  2A3  -  A2  +  8A) 

and  —  2A3  +  A2  —  8 A  is  not  the  square  of  some  polynomial  in  C[A],  This  shows 
that  P(A,  M)  is  irreducible  as  a  quadratic  polynomial  in  M  over  C[A],  so  in 
particular  P(A,  M)  is  irreducible  in  C[A,  M].  (A  similar  argument  shows 
that  the  polynomial  Q( A,  M )  =  2A 3M2  +  A2(12M  +  4)  +  8A  +  4  of  Case  2  is 
also  irreducible.) 

Assume  that  m  is  an  integer  such  that  P( A,  m)  is  still  irreducible  as  a 
polynomial  of  degree  3  in  A  with  integer  coefficients.  Let  9rn  be  any  root  of 
P(A,  m)  and  let  Am  be  the  discriminant  of  the  cubic  field  Q[#m],  Corollary 
3.12  in  [3]  shows  that  the  number  of  integer  solutions  (bb,  A)  of  equation  (8) 

is 

«  |Am|°'201.  (9) 

Notice  that  Am  is  a  divisor  of  the  discriminant  of  the  polynomial  A3  +  4A2  — 
(12m  —  4)A  +  4m2  (obtained  by  rewriting  the  equation  P(A,m)  =  0  as  a 
monic  polynomial  equation  in  2/A)  and  this  last  discriminant  is  — 16 (27m4  — 
216m3  +  280m2  —  48m).  (In  Case  2,  Q( A,  m)  produces  the  polynomial  A3  + 
4A2  +  (12m  +  4)A  +  4m2,  whose  discriminant  is  — 16(27m4  + 216m3  +  280m2  + 
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48m).)  In  both  cases,  these  expressions  are  never  zero  if  m  is  a  non- zero 
integer,  which  is  an  observation  that  will  be  useful  later.  For  the  moment, 
we  simply  record  that  estimate  (9)  and  the  above  calculation  show  that  the 
number  of  integer  solutions  to  equation  (8)  is  <C  |m|'804.  Clearly,  every  integer 
solution  (IF,  A)  of  equation  (8)  determines  p  uniquely.  This  shows  that  the 
number  of  primes  p  that  can  arise  in  this  way  when  P( A,  m)  is  irreducible 
and  1  <  Iml  <  Y  is 


« 


E  ■ 

1<M<Y 


.804 


<  Y 


1.804 


•c 


Y 1-804 
J\902  ' 


Assume  finally  that  1  <  |m|  <  Y  is  such  that  P(A,  m)  is  reducible  as  a 
polynomial  in  A  with  integer  coefficients.  By  Hilbert’s  Irreducibility  Theo¬ 
rem,  the  number  of  such  values  for  m  is  <C  Y 1//2  log  Y  <C  F1/2  log  A”  (see,  for 
example,  [2],  or  Theorem  1,  Section  13.1  in  [7]).  Multiplying  both  sides  of 
equation  (8)  by  4m4,  we  get 

U2  =  F3  -  (12m  -  4)F2  +  16 m2V  +  16m4, 

where  U  =  2m2W  and  V  =  2A m2  (in  Case  2,  we  get  the  equation  U 2  = 
F3  +  (12m  +  4)V 2  +  16m2F  +  16m4  with  similarly  defined  U,  V,  as  in  Case 
1).  Write 

A3  —  (12  m  —  4)  A2  +  16  m2  A  +  16  m4  =  (A  +  a)  (A2  +  bA  +  c), 


where  a,  b  and  c  are  integers.  Then  there  exist  two  squarefree  integers  d\ 
and  d2  dividing  the  discriminant  of  the  above  polynomial  (which  is  16m4  Am), 
such  that 

V  +  a  =  d^U2  and  V2  +  bV  +  c  =  d2U22. 

The  second  equation  can  be  rewritten  as 

V2  +  A  =  d2Ul  (10) 

where  V\  —  2V  +  b,  U3  =  2 U2  and  A  =  4c  —  b2.  Note  that  A  /  0,  since 
if  A  =  0,  then  P(A,  m)  has  a  double  root,  and,  as  we  have  seen,  this  is 
not  possible  if  m  ^  0  is  an  integer.  Note  that  ac  =  16m4  and  ab  +  c  = 
16m2,  therefore  |c|  <  16m4  and  |6|  <  16m2  +  16m4  <  32m4,  which  shows 
that  |A|  -C  m8.  It  is  now  easy  to  check  that  max{|Fi|,  |C/3|}  -C  A"8.  The 
arguments  from  the  first  part  of  our  proof  (the  case  when  A  <  T)  show 
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that  for  fixed  values  of  A  and  d2,  the  number  of  integer  solutions  (Vi,U3) 
to  equation  (10)  is  -C  r(|A|)logA  <C  as  A"  — >  oo.  Since  d2  can 

be  chosen  in  at  most  r(|A.m|)  <  X0^  ways  for  a  fixed  m,  there  are  only 
0(Y^2  log  A)  possibilities  for  m,  and  each  quadruple  (m,  d2,  V\,  U3)  arising 
in  this  way  determines  p  uniquely,  we  get  that  the  number  of  possibilities  for 
p  is  at  most  Ylt2X°^  <  X 2/3  whenever  A"  is  sufficiently  large.  (The  same 
argument  holds  for  Case  2.)  Putting  everything  together,  we  get  that  the 
number  of  primes  p  G  (A/2,  X)  such  that  the  Diophantine  equation  (1)  can 
have  a  non-trivial  proper  solution  is 

Y 1-804 

«T(logA')2+  y^+XW. 

Choosing  T  =  x1’804/1-902,  we  get  that  the  number  of  such  possibilities  for  p  is 
<  A0'95  for  large  values  of  X.  Since  there  are  >  (0.5  +  o(l))Ay  log  A"  primes 
p  G  (A/2,  A)  as  X  — >  oo,  we  deduce  that  for  most  primes  p  in  the  above 
interval,  equation  (1)  has  no  non-trivial  reduced  solutions  when  N  =  2 p. 
This  completes  the  proof  of  the  theorem. 

3  Related  elliptic  curves 

In  the  construction  of  solutions,  Stroeker  obtained  elliptic  curves  of  the  fol¬ 
lowing  equations  (see  formulas  (14)  and  (14’)  in  [6]) 

^2  =  4(A  -  1) (v  +  l)3  +  (A  -  3v  -  2)2,  (11) 

and 

w2  =  4(A  +  l)(w  —  l)3  +  (A  —  3w  +  2)2.  (12) 

If  u  =  (A  +  3)  / 4,  then 

{viz)  —  ((A  —  3)/4,  A(A  +  l)/4)  (13) 

is  a  point  on  the  elliptic  curve  dehned  by  (11).  In  the  same  way, 

(u,  w)  —  ((A  +  3)/4,  A(A  —  l)/4)  (14) 

is  a  point  on  the  elliptic  curve  defined  by  (12). 

Theorem  2.  The  above  points  ( v ,  z)  and  ( u ,  w )  are  not  torsion  points  except 
if  A  =  3,  9.  Moreover,  if  A  =  3,  both  points  are  of  order  9  and  if  A  =  9, 
they  are  of  order  6. 
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Proof.  Using  Magma  [1],  for  2  <  N  <  100,  we  checked  that  the  points  (13) 
and  (14)  are  not  torsion  points  on  the  curves  (11)  and  (12),  respectively, 
except  for  N  =  3,9.  If  N  =  3,  both  are  torsion  points  of  order  9,  while 
when  N  =  9,  both  are  torsion  points  of  order  6.  To  achieve  this,  multiply 
both  members  of  equation  (11)  by  4(1V  —  l)2  and  put  X  =  4(N  —  l)v  and 
Y  =  4(N  —  1  )z  to  obtain 

Y2  =  X3  +  (127V  -  3)X2  +  24N(N  -  1)X  +  (4N(N  -  l))2. 

Now  we  set  y  =  Y  and  x  —  X  +  41V  —  1  and  arrive  at 

y2  =  x3  +  (— 241V2  -  3)x  +  161V4  +  401V2  -  2.  (15) 

If  u  —  (N  +  3)/4,  then  v  —  (IV  —  3)/4  and  z  =  N(N  +  1) / 4.  One  can 
deduce  that  x  =  N(N  —  1)  and  y  =  5 N  —  4.  Therefore  P  =  (x,  y)  =  (IV2  + 
2,  N(N2  —  1))  is  the  point  on  (15)  corresponding  to  the  initial  point  (v,  z )  = 
((IV  —  3) / 4,  N(N—  1) / 4)  on  (11)  via  the  above  birational  transformation.  We 
follow  a  similar  procedure  for  equation  (12)  and  we  obtain  the  same  equation 
(15)  with  the  same  corresponding  point  P  =  ( x ,  y)  =  (IV2  +  2,  N (IV2  —  1))  to 
the  point  (u,w)  when  u  =  (N  +  3)/4  on  (12).  So,  we  deal  only  with  equation 
(15). 

Let  kP  =  (. xk,yk )  be  the  sum  of  P  =  (x,y)  =  (IV2  +  2 ,N(N2  —  1)) 
with  itself  k  times  in  the  Mordcll-Weil  group  of  (15).  By  Mazur’s  Theorem 
(see  [5]),  if  P  is  a  torsion  point  then  kP  =  O  for  some  k  G  {1, . . . ,  12}. 
So,  we  computed  ( Xk,yk )  for  all  1  <  k  <  12.  In  fact,  consider  the  associated 
projective  points  (Xf..  :  Y*.  :  Zjf)  with  Xk,  W,  Zk  G  Z[!V],  We  solved  separately 
each  one  of  the  polynomial  equations  Yk(N)  =  0,  Zk(N)  —  0  for  2  <  k  <  12. 
Here  are  the  results: 

•  If  Yk(N)Zk(N )  =  0  for  some  positive  integer  N  and  1  <  k  <  12,  then 
k  G  {3,6,9,12}. 

•  If  Y:i(N)  =  0,  then  N  =  9; 

•  If  Z6(N )  =  0,  then  N  =  9.  If  Y6(N)  =  0,  then  N  =  9; 

•  If  Z9(N)  =  0,  then  N  =  3; 

•  If  Z12(N)  =  0,  then  N  =  9. 
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